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Abstract. The electron transport properties of ZnO nano- wires coupled by two 
aluminium electrodes were studied by ab initio method based on non-equilibrium 
Green's function approach and density functional theory. A clearly rectifying current- 
voltage characteristics was observed. It was found that the contact interfaces between 
Al-0 and Al-Zn play important roles in the charge transport at low bias voltage and 
give very asymmetric I-V characteristics. When the bias voltage increases, the negative 
differential resistance occurs at negative bias voltage. The charge accumulation was 
calculated and its behavior was found to be well correlated with the I-V characteristics. 
We have also calculated the electrochemical capacitance which exhibits three plateaus 
at different bias voltages which may have potential device application. 
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1. Introduction 

The investigations on ZnO material have been continued for many decades due to its 
wide range of apphcations [1]. Due to its large direct band gap {Eg ~ 3.3 eV at 300 
K) and large exciton binding energy ~ 60 meV (room-temperature 300 K is about 
~ 25 meV), it is an important material for manufacturing ultraviolet devices [2] and 
high temperature luminescent devices [31 H]. Its easy growth with high purity makes 
it a low cost material as high quality substrate [5l E] and its easily etched property 
makes it possible to be made as small-size devices. It can be used as space devices [7] 
and transparent thin-film transistors [8], [9] because of its insensitivity to high energy 
radiation and visible light. It can also be a candidate as a spintronic material because 
Mn-doped ZnO, as well as Fe-, Co-, or Ni-alloyed ZnO, has high Curie temperature 
ferromagnetism [10]. There are many other properties and applications for ZnO bulk 
material and ZnO related materials that make ZnO a unique material for emerging 
nano-technology [1]. 

Recently, several configurations of ZnO nanostructures such as nanowires 
nanotubes [121 [131 [H], nanobelts [151 HH], nanosheets [171 HH] and nanorods [121 fT4] . 
have attracted much attention due to their potential applications and easy fabrication 
[T9] . Such ZnO nanostructures have novel electrical [201 1211 [22] and optoelectrical 
[231 [21] properties as well as mechanical properties [19]. They can be used to make 
sensitive chemical nanosensors [251 [26], solar cell [27], light-emitting devices [291 [28] , 
nanogenerators (converting nano-scale mechanical energy into electricity) and nano- 
piezotronic devices [19]. First principles calculations for the geometry and electronic 
structure of small ZnO nanowire and nanobelt were also carried out by many authors 
[301 [311 [32] . More recently, the carrier transport of ZnO coupled by metal layers such as 
Au electrode and Mg electrode were studied by Kamiya et al [33] • In their calculation, it 
was found that the Au-ZnO-Au interface shows Schottky contact behavior while the Mg- 
ZnO-Mg interface shows ohmic contact. The contact behavior has crucial influence on 
the device property and is very sensitive to the fabrication processes [34] . Hence it is very 
useful to study the contact behaviors between ZnO and other materials. In addition, 
it is very useful to study other transport properties such as charge polarization and 
electrochemical capacitance. In this paper, we concentrated on the transport behaviors 
and charge distribution and charge accumulation of Al-ZnO-Al molecular structure by 
the density functional theory combined with nonequilibrium Green's function method. 
Different from Au-ZnO-Au and Mg-ZnO-Mg nanostructures, we found that for Al-ZnO- 
Al molecular junction the I-V curve is highly asymmetric. In addition, a negative 
differential resistance appears at negative bias voltage. Furthermore, we have calculated 
the electrochemical capacitance of this device and found that depends on bias 
voltage in a unique way: there exists three regions where is constant. At small 
bias, our analysis shows that electrochemical capacitance is dominated by the classical 
capacitance whereas at large bias voltage quantum effect becomes important. 
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2. Method 



Here we present our first principles analysis of electronic transport through a Al-ZnO- 
Al nano-device. The calculation package we used is NADCAL [35] (the developed code 
of McDcal [36]), which is an implementation of real space density functional theory 
(DFT) within the Keldysh non-equilibrium Green's function formalism (NEGF). The 
DFT method provides the system electronic structure and NEGF includes information 
of the non-equilibrium quantum statistics. The system we studied was a typical two- 
probe open system which was divided into a scattering region and two semi-infinite 
electrodes. The electronic structure of the central region was calculated self-consistently 
and the couplings between the central region and the two semi-infinite electrodes were 
included by the self-energy. Enough buffer layers for Al atomic electrode were included 
in the cental region so that the boundary conditions of the scattering region could be 
determined through the corresponding infinite left and right electrodes. We refer to 
references [36l ETJ [38] for detailed description of NEGF-DFT technique. 

Our central scattering region consisted of 6-layers of (0001) ZnO with width 3.752 
A and three buffer layers of each aluminium electrodes (see figured]). Each ZnO layer 




Figure 1. Schematic plot of Al-ZnO-Al (Al: gray; Zn: white; O: black). The direction 
from left to right is considered as z-axis. 
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Figure 2. Contact distances between aluminium electrode and ZnO determined by 
finding the lowest energy, (a) Total energy versus the distance between Al-0 contact; 
(b) Total energy versus the distance between Zn-Al contact. 
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had 3 oxygen atoms and 3 zinc atoms. Basis set was composed of s,p,d real space 
linear combination atomic orbitals (LCAO). 3(i^°4s^ was chosen as valence electrons 
for zinc and 2s^2p^ for oxygen. Exchange correlation was treated using LDA_PZ |39j . 
When solving the Poisson equation in real space one grid point occupied a volume of 
4.583 X 10"'^ A^. In the self-consistent calculation, the tolerance was set to be less than 
10~^ Hartree. 

For Al-ZnO-Al structure in figure [H the left electrode was coupled directly to 
oxygen atoms (Al-0 contact) and the right electrode was coupled to zinc atoms (Al- 
Zn contact). Instead of relaxing ZnO and aluminium electrodes, we determined the 
distances between ZnO and electrodes by finding the minimum of total energy when 
changing the distance of Al-0 and Zn-Al respectively. According to figure [2] we set the 
distance between aluminum electrode and oxygen side of ZnO to be 1.95 A while the 
distance of Zn-Al to be 2.75 A. 

For the two-probe NEGF-I-DFT formalism, the electronic current can be obtained 
by Landauer-Biittiker formula 

/ = ^| dEi:T[f{E,VL.VR)]UL{E,iiL) - ^R{E,^^n)), (1) 

where T{E, Vl, Vr) is the transmission matrix and is obtained by Green's function from 
the following formula 

f{E,VL,VR) = TLG''TRG''. (2) 

Il/r = [1 + exp((£' — Hljr) /kT)]~^ is the Fermi distribution function of left/right 
electron reservoir connected to the left/right electrode at infinite z = — oo/oo. ^l/r = 
fJ'i/r + QVL/R is the electrochemical potential with /i//,, the chemical potential of left /right 
lead and Vl/r the apphed potential on the left /right lead. 

3. Results and Discussions 

In figure |3l we plotted the current-voltage (I-V) curve of Al-ZnO-Al. A rectifying I-V 
characteristics is clearly observed. The current increases slowly with increasing of the 
positive bias voltage Vh but decreases quickly when the bias voltage reverses its sign. 
Different behaviors of I-V characteristics when the bias voltage changes from positive 
to negative can be understood as follows. At zero temperature and in hartree unit, 
equation (|T]) becomes (We set Vr = in our calculation.) 

/■o dE 

1 = 2 ^TrT(E,H). (3) 

We find that when the positive bias voltage is applied, the current is contributed 
from the energy integral interval E = (— |Vfe|,0). But for negative voltage H < 0, 
the transmission coefficient will be integrated on the positive interval E = (0, |Vfe|)- 
Although the transmission coefficient is a function of the bias voltage, we first analyze 
the I-V curve near Vh = 0. For this purpose, we plotted the transmission coefficient 
and the density of states (DOS) versus the Fermi energy at H = in figure [3l In 
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our calculation we find that three kinds of states give major contributions to the I-V 
curve. They are zinc 4s-orbital, oxygen ^-component of 2p-orbital and oxygen 2s-orbital. 
Among them, contribution due to oxygen 2s-orbital is the smallest so we focus on the 
other two orbits. It is clear that when E > Q the DOS of zinc s-orbital is much larger 
than that of oxygen p^-orbital and when < both zinc s-orbital and oxygen p^-orbital 
contribute near E = Q but the contribution is less than that of the zinc s-state at > 0. 
So we conclude that when negative bias voltage is added electrons pass device mainly 
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Figure 3. (a) Current as a function of bias voltage Vb = Vl — Vr; (b) Transmission 
coefficient versus energy at bias voltage VJ, = 0; (c) DOS versus energy due to zinc 4s 
state (Zn:s solid line), oxygen 2p state in z-direction (0:pz dashed line) and oxygen 2s 
state (0:s dotted line) 
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Figure 4. (a) Local density of states of zinc 4s-state on each layer; (b) Local density 
of states of oxygen pz-state on each layer. 



Electron transport through Al-ZnO-Al: an ah initio calculation 



6 



from zinc s-orbital but when the positive bias voltage is added electrons pass device 
mainly from both zinc s-orbital and oxygen p^-orbital. To provide further evidence, 
we analyze the local density of states (LDOS) according to the layers of ZnO in the 
scattering region (see figure H]). Note that we label the layers of oxygen and zinc from 
left to right. One can find that each layer has very different LDOS. For zinc atoms, 
the rightmost layer (the layer 6, i.e. the layer coupled directly to aluminium lead) gives 
the largest LDOS around the Fermi level. (The Fermi level is shifted to = in 
our calculation.) When the energy increases the LDOS of middle layers (layers 2 to 4) 
gradually increase and reach their maxima at about E = 0.83 eV. But the LDOS of 
all zinc layers become small when E < —0.25 eV. For oxygen atoms, similar result is 
found if we recall that the first oxygen layer is connected directly to the left aluminium 
electrode. The LDOS of oxygen give the largest value at the first oxygen layer around 
the Fermi level E = 0. When E becomes quite small and closes to £' = — LO eV, 
the LDOS of the middle layers (layer 3 and layer 4) dominate and when E > the 
LDOS of oxygen are small compared with the those of £' < 0. In another word, when 
small negative bias voltage is applied the contact interface between Al-Zn dominates 
I-V behavior. But when the positive bias voltage is applied, both Al-Zn interface and 
Al-0 interface will dominate the charge transport. The rectifying behavior is mainly 
controlled by the contact interface characteristics due to the different coupling atoms 
Al-Zn and Al-0 when bias voltage is small. 

Actually the transmission coefficient T = TT[T{E,Vb)] is a function of the bias 
voltage Vb = Vl — Vr at zero temperature so the situation becomes more complicated 
at large VJ,. In figure O the transmission coefficient as a function of the energy E and 
the bias voltage H was depicted. The I-V curve for positive bias voltage corresponds 
to the triangular area A between the vertical line E = and the slope line Vb = —E 
for bias voltage H > and energy E < (see the dashed black lines in figure [5]) and 
the I-V curve for negative bias voltage corresponds to the triangular area B between 
the above two lines for Vf, < and E > 0. We find that in the right low triangular 
area B the transmission coefficient is much large within the interval Vb G (—0.4, —0.2) 
V and reach maximum around Vb ~ —0.35 V (see the C region in figure [5]). That's why 
the absolute current value decreases when the absolute value of bias voltage increases 
around Vb G (—0.5, —0.35) V and the negative differential resistance appears in that 
voltage region (see I-V curve on the interval (—0.5, —0.35) V in figure [3](a)). 

For the quantum devices, due to finite density of states the charge accumulation 
in the scattering region differs from the charge accumulation in the classical metal 
capacitance and one should use the concept of electrochemical capacitance = 
AQ/A/i, where A/i is the electrochemical potential difference between the two electron 
reservoirs jlQl HH |12] and the AQ is the charge accumulation. For our two-probe 
problem, both electrodes are aluminium so that fii = fir- Hence, at zero temperature, we 
have Afi = Vb = V^ — Vr. When the bias voltage is applied, the electrochemical potential 
across two probes is changed by A/i = Vb from equilibrium, charges will be injected 
into the scattering region. The induced charges will occur due to the electron-electron 
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Figure 5. Transmission coefficient versus bias voltage VJ, and energy E. The dashed 
black lines are only used as symbol lines in order to depict the region A, B and C. 



interaction. In DFT, the injected and induced charges are automatically considered 
by solving the Kohn-Sham equation and the Poisson equation self-consistently and the 
Gauge invariant is also satisfied automatically |13] . Such a charge accumulation is closely 
related to the I-V characteristics and can only be determined by solving scattering 
problem. For a rough estimate, we divide the system into two parts from the middle 
of scattering region called left region and right region. We then calculate the charge 
accumulation A(5l(H) = — Ql(0) and A(5_R(Vb) = Q_R,(Vfe)— Qi?(0) and plot them 

in figure El The electrochemical differential capacitance = dQ/dVb was also plotted 
in figure [61 Here, Q is the fitted curve of AQjiiVb). Following observations are in order. 
First of all, the electrochemical capacitance depends on the voltage although it is a piece- 
wise constant in three regions labeled as (a), (b) and (c). Secondly, the asymmetry 
of electrochemical capacitance is clearly seen. In region (a) which is asymmetric 
with respect to the bias voltage H = (—0.27, 0.01) V the electrochemical capacitance 
assumes a smaller constant value = 0.003 aF. While in region (b) and (c), are 
approximately four times larger. The value of the electrochemical capacitance also 
exhibits asymmetry between region (b) and (c). Finally, the behavior of electrochemical 
capacitance in region (a) can be understood qualitatively as follow. Since in region (a) 
the current is small, we can use the following formula for electrochemical capacitance 
(where the current is assumed to be zero) [H] 



Co cWl eWn' ' ' 

where Cq is the classical capacitance which can be roughly estimated using the parallel 
plate capacitor Co = eoA/d ~ 0.002 aF where A is the cross-section area and d is the 
separation of two plates. Here and are the DOS of the left and right regions, 
respectively. From figure HI we see that Dl and Dr are much larger than Co/e^ and 
hence can be neglected in equation (jl]). Therefore, we conclude that in region (a), the 
electrochemical capacitance is dominated by the classical capacitance. In region (b) and 
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Figure 6. (A)Charge accumulations versus bias voltage on left scattering region 
{AQl- triangles) and right scattering region (AQi?,: squares). The black solid line 
Q{Vb) is the fitted curve of AQu; (B) Electrochemical differential capacitance dQ/dVb 
versus bias voltage. 

(c), however, quantum effect must be considered. At the moment, there is no analytic 
formula to account for the electrochemical capacitance at the nonlinear bias voltage. 

In summary, we have analyzed the electron transport through the Al-ZnO-Al atomic 
junction and clear rectifying characteristic was observed. It was found that the contact 
interfaces between aluminium electrodes and central device ZnO play essential roles 
in the current-voltage behaviors when the bias voltage is small. When the negative 
bias voltage is applied, the major contribution to the electron transport comes from 
Zn s-states. But when the positive bias voltage is applied, both Zn s-state and O p^- 
state are responsible in the charge transport. The negative differential resistance was 
found in our device. Our results showed that for such a quantum open system the 
charge accumulations in the scattering region are closely correlated with the applied 
bias voltage and I-V characteristics. The electrochemical capacitance was found to have 
three plateaus. At small bias, the electrochemical capacitance is dominated by the 
classical capacitance. At large bias, however, quantum effect becomes important. Our 
findings may have potential device application for future nano-technology. 
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